Extended Standard Map with Spatio- Temporal Asymmetry 
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We analyze the transport properties of a set of symmetry-breaking extensions of the Chirikov- 
Taylor Map. The spatial and temporal asymmetries resuh in the loss of periodicity in momentum 
direction in the phase space dynamics, enabling the asymmetric diffusion which is the origin of the 
unidirectional motion. The simplicity of the model makes the calculation of the global dynamical 
properties of the system feasible both in phase space and in controlling-parameter space. We present 
the results of numerical experiments which show the intricate dependence of the asymmetric diffusion 
to the controlling parameters. 
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I. INTRODUCTION 

The motion of a particle in temporally and spatially 
periodic potential has attracted renewed attention as a 
model of molecular motors Q . In order to make it a mo- 
tor, that generates unidirectional motion out of periodic 
motion, the potential has to be a "ratchet" , i. e. its reflec- 
tion symmetry in both spatial and temporal direction has 
to be broken, a fact generally known as a Curie's theorem. 
Traditionally, the friction has supplied the source of tem- 
poral symmetry breaking ^, |^, ^ . There are several re- 
cent works 1^, |7| in which the frictionless "inertial" regime 
is studied where the time-asymmetric perturbation is the 
driving force of the unidirectional motion. Through these 
studies, it has become evident that the characeristics of 
the phase space plays a dicisive role. The exploration 
of the full phase space of a dynamical system, however, 
requires considerable computational resources. This nat- 
urally leads to models of ratchet dynamics in which the 
differential equation governing the time evolution is re- 
ducible to the discrete difference equations, namely the 
map . One would then like to have a model of a system 
with ratchet dynamics, which is free from ad hoc assump- 
tions and parameters that tend to obscure the underlying 
simple physics. 

In this article, we propose a simple model of a particle 
motion under temporally and spatially periodic potential 
whose key feature is in that the obtained descrete map 
is a natural extention of the Chirikov- Taylor map, also 
known as the standard map Q whose dynamical proper- 
ties have been thouroughly analyzed. The unidirectional 
motion of the system is understood in the context of the 
deterministic diffusion foimd in the standard map. Also, 
as a practical benefit of the discrete evolution of the map 
dynamics, the computational burden is substantially re- 
duced, and it becomes possible to study, for the first time, 
the transportation properties of a system in a broad re- 
gion in the controlling parameter space. A very intricate 
dependence of the unidirectional motion on the control- 



ling parameter is revealed though a numerical investiga- 
tion. 



II. EXTENDED STANDARD MAP: 
FRICTIONLESS CASE 

We consider the motion of a classical particle described 
by a time-periodic Hamiltonian with instantaneous po- 
tential action, or (5-function kick, which has space peri- 
odic part and alternating-current type driving motion; 

2 ^ 

H=Y+ [u{x)+xAsn]5{t-nT{l + eSn}), (1) 

n— — oo 

where, for now, the trapping potential u{x) is set to be 
symmetric and given by 

u{x) — —Kcosx, (2) 

and Sn = (—1)" is the parity of the iteration number n. 
The evolution of the system is described by 

X = p (3) 

OO 

P = - W{x)+Asn]5{t-nT{l + es„}), (4) 

n— — oo 

One can integrate the motion in between the ^-function 
kicks and obtain the map which describe the evolution 
of position Xn and momentum p„ at discrete time = 
nT{l + es„) - 0; 

Xn+i a;„ +p„+i(l-|-es„) (5) 

Pn+l = Pn - K sinXn - ASn- (6) 

This set of equations can be regarded as a three- 
parameter extension of the standard map The pa- 
rameter K, which is present in original standard map 
is the strength of the periodic potential that traps the 
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particles in each periodic location. The parameter A is 
the amplitude of alternating swing motion, and e is the 
measure of its time asymmetry. As with the case of the 
standard map, the dynamical properties of the system 
is studied through the analysis of phase space The 
transport property of standard map, A = e = 0, is char- 
acterized in terms of phase space average of functions of 
position and momentum x and p. It is well known that 
the direct averages {x) and {p) are zeroas a direct result of 
the fact that the standard map is both time-reversal and 
space-reflectional symmetric. However the mean square 
(p^) can be non zero when one has the chaotic motion un- 
bounded by KAM tori in p-direction. This happens with 
sufficiently large value of K. One then has deterministic 
diffusion through chaotic motion which is characterized 
by the random- walk behavior, (a;^) oc t. In the extended 
standard map (^)-@, a unidirectional motion can be re- 
alized if one has by the non-zero value for the (p) itself 
in the diffusive motion. 

The reflection symmetries play an essential role in the 
existence (or non-existence) of the unidirectional motion 
1^. One can rewrite (||)-@ as 

- Xn+i = -a;„ -p„+i(l + £s„) (7) 

-Pn+l = -Pn - K sin {-Xn) + ASn. (8) 

These become identical in form to (||)-@ themselves in 
case oi A = 0. Therefore, in this case, one has 

(x„) = (-x„) = 0, (9) 

{Pn) = {-Pn)=0 (forA = 0), 

which means that there cannot be any directed transport. 

We now define x„ = x^ and p„ = Pn+i, obtain a re- 
verse map 

Xn-l = i„ -p„-i(l - £s„) (10) 

-Pn-l = -pn ~ u{Xn) ~ ASn- (H) 

Therefore, if one sets s — 0, the map for {a;„,p„} and 
{x-m —p-n} are identical, and one has 

i-p) = i-p) = {p)=0 (for e = 0) (12) 

Therefore, one needs both A and e ^ to have 
directed motion. It is important to note that the map 
(^"® ™ f^*^* differently for odd and ven n because 
of the existence of Sn term. It is advantageous to consider 
the phase spaces for odd n and even n separately. The 
reflection symmetry between (H)-(||) and (0)-(||) should 
then be interpreted as one between two phase space pro- 
files each representing odd n and even n. 

For the standard map, A = e = 0, as is well known, 
the map is 27r-periodic in both x and p directions. The 
periodicity for x direction is kept intact for non-zero A 
and e: when one has a map (x„,p„) — > {xn+i,Pn+i), one 
also has 

{Xn + 2Tr,Pn) (Xn+l +2Tr,Pn+l). (13) 
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FIG. 1: Phase space profiles of the map Eqs. (5)-(6) for 
even numbers for n. The strength of periodic potential is 
set to be A' = 0.6. The temporal asymmetry parameter e 
and amplitude of alternating swing A are varied. Left (right) 
column is for odd (even) n. 

On the other hand, original 2tt periodicty for p - direction 
is lost when non-zero e is present. If, however, e is a 
rational number, M/N where M and N are mutually 
incomensurate positive integers, a new 2iV7r-periodicty 
for p-direction emerges, because in this case, one has 

{Xn,Pn + 2N7T) ^ (14) 
{N ±M)TT,Pn+l+2NTT) 

(for e = M/N;N > M > 0). 
The change of the p-periodicity has a direct impact on 
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FIG. 2: Phase space orbits of the map (5)-(6) for e = (top) 
and e — 0.2 (bottom) started from a set of random initial 
conditions confined in the region \p\ < vr/lO and |a;| < n/2. 
The strength parameters are set to be A = 1.2 and K — 0.6. 
For the case of e = 0, one has (p)^^^ + (p)even = 0- 



the transportation properties of the extended standard 
map. The 2iV7r-periodicity in p-direction for a rational 
number e — M/N (where M and N are two mutually 
incommensurate natural numbers) means that the KAM- 
tori (when they exit) appear at least once in 2Nt: period. 
Thus the chaotic diffusion is bounded within the region in 
between two adjacent tori. This fact gives an estimate for 
the average momentum {p) within a region that include 
p = in the form 



(p) < Ntt (for £ = M/N; N > M > 0). 



(15) 



Note that only the upperbound is obtained since the ex- 
act location of the two adjacent KAM tori near p — 
cannot be determined a priori. This is evident when one 
considers s = case, in which the map is 27r-periodic in 
p-direction, and a similar estimate to above gives (p) < 
TT, while, as we have lernt in (|l2|), it is zero because KAM 
tori is placed symmetrically around p = 0. 

The situation is illustrated in Fig.l where we draw 
phase space profiles for the frictionless extende standard 
map with various e and A. The strength parameter is 
chosen to he K — 0.6 in these examples. The effect of 
the breaking of phase space symmetries on the trans- 
portation properties of the system becomes evident by 
the inspection of Fig. 2, where we plot the portion of 
the phase space profile that is reached from a set of ini- 
tial conditions clustered around in the vicinity ot p — 0. 
With the onset of non-zero A and non-zero e terms, the 
average momenta in the chaotic region around p = is 
shifted from zero, and the directed tranport is realized. 
This is essentially in accordance with the conclusion of 
Flach et al. in which they identify this phenomenon 
asymmetrical Levy flight. We are now ready to ready to 
look at the directed transport not just in one parameter 




FIG. 3: A measure of directed motion (p) as the function of 
asymmetry parameter e and the strength of swing motion A 
for the extended standard map with asymmetric kick (5)-(6) 



set, but in a global region in the parameter space. In 
Fig. 3, we plot the average momentum {p) of the chaotic 
band that include p = as the function of swing param- 
eter A and time asymmetry s. The coupling strength 
K is taken to he K = 0.5 (top), K — 0.6 (middle) and 
K = 0.7 (bottom). We start from a ensemble of 800 ran- 
dom initial configuration satisfying \p\ < tt/IO and |a;| < 
7r/2. We symmetrize the ensemble in terms of p— and x— 
to guarantee (p) = exactly at the start. Each config- 
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uration is evolved with (5)-(5) for 12,000 times, and the 
momentum p is averaged over that period and over the 
ensemble to obtain (p) for a given A and e. Calculation 
on the grid points of 48 x 42 in the space of {A, e) is used 
for each of the three figures. 

From these figures, it appears that the direction of the 
asymmetric transport (positivity or negativity) is con- 
trollable, having step-like structure of increasing (p) with 
decreasing e, which is in broad accordance to our ar- 
guments leading to ([1^). An intriguing feature is the 
direction-reversal of (p) with parameter variation, whose 
occurence is not at all obvious in the outset. Note also 
that the speed of unidirectional transport can be rather 
large: {p) is almost 27r in the peak region, which means 
that the particle moves one spatial period in one tempo- 
ral period. 

III. RATCHET MAP WITH FRICTION 

Our model can be extended to incorporate the ratchet 
dynamics with damping through the addition of friction. 
In place of eqs. (^)-@, we set 

X ^ p (16) 

oo 

P ^ - ["As(a;)+^s„] (17) 

n— — cxj 

y.5{t - nT{l + es„}) - ^p 

where we now add spatially asymmetric term to the trap- 
ping potential as 

u'j^g{x) — K{smXn + 2/xsin(2a;„ -f- 5)}. (18) 

The map obtained from these is written as 

xn+i = x„+p„+i(e^(i+"")-l)/7 (19) 
Pn+i = b„ - u:is(a:„) - As„] e-^(i+"") (20) 

There are now six parameters in the model: K, A, fi, S, 
7, and e. We reduce this number by setting e = 0. Since 
the temporal asymmetry is already broken by 7 alone, 
we do not expect to loose out much. We now obtain 

Xn+l = Xn+Pn+l{e^ (21) 
Pn+1 = [Pn - u'j^g{Xn) - ASn] e"^ (22) 

As in the frictionless case, we consider the reflection 
symmetry of the map (pT|)-(p^). With with the diffinition 
in = Xn and pn = Pn+i, ouc obtaius a reverse map 

i„„i = x„ -p„_i(7 - l)/7 (23) 
-Pn-i = -_p„e'^ + [-u^5(£„) - As„]e"'^ (24) 

In the elastic limit 7 — > 0, the map for {x„,p„} and 
{x-m —p-n} are identical, and one has 

{~p) = {~p) = (p) = (for 7^0) (25) 
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FIG. 4: Phase space orbit of map (21)-(22) after 10000 iter- 
ation showing the starnge attractor for A = 1.8 (top figures), 
and the periodic attractor A — 2.0 (bottom figures). Other 
parameters are set to he fi = 0.2 and S = 0.3. 

irrespective to the symmetry property of uas- This re- 
sult, puzzling in first inspection, is of course in accor- 
dance with Curie's theorem. If, on the other hand, we 
have A — 0, the motion is ovcrdamped, and in sufficient 
steps, it will come to stand still, namely, 

x^xo {u'asM^O), p^O (forA-0) (26) 

With friction, after sufficiently large steps n, all trajecto- 
ries fall into the atractors. They can be either stable at- 
tractors with finite period, or the strange attractors with 
fractal dimension depending on the parameter value [ p4^ . 
Two examples from each case are shown as the top and 
bottom figures in Figs. 4. As is well known, these two 
cases alternate with minute variation of parameters. If 
one has n = 5 = 0, the shape of attractor respects the 
spatial reflection symmetry, when both even and odd n 
are considered. Therefore, one can only have {p) in this 
case. On the other hand, as long as /i 7^ and S ^ 0, one 
has no spatial mirror symmetry in the phase space, and 
one should have (p) 7^ in general. It is expected that 
the value of (p) can vary wildly as fuctions of parameter 
because of the afore-mntioned alternation of periodic and 
strange attractors with minute change of he parameter. 

The transport property of the extended standard map 
with friction can be overviewd in Fig. 5, in which we 
show the three-dimensional plot of (p) as functions of two 
controlling parameters A and 7. The spatial asymmetry 
parameters are fixed to be /i — 0.2 and S — 0.3. The 
coupling strength K is varied from iir=:1.0toiir=1.2 
from the top to bottom. As before, we start from a p— 
and x— symmetrized ensemble of 1600 initial configura- 
tion satisfying |p| < tt/IO and \x\ < 7r/2. Each config- 
uration is evolved with (pT|)-(^) for 40,000 times, and 
the momentum p is averaged over that period and over 
the ensemble to obtain (p) for a given A and e. The grid 
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FIG. 5: The average (p) as the function of 7 and A for the 
extended standard map with damping (21)-(22) The param- 
eters controling the spatial asymmetry are set to be = 0.2 
and 5 — 0.3. 



points of 64 x 64 in the parameter space of {A, 7) in the 
range A e [0.8,2.8] and 7 G [0.02,0.10] are used for each 
of the three figures. 

One can see from these figures, that the overdamped 
motion in ratchet potential shows directed transport, but 
is less controllable than the inertial counterpart, as ex- 
pected from our preceeding arguments. At the moment, 
we do not have any clue as to the physical origin of the 
intriguing structure found in Fig. 5. However, we feel 
that the simplicity of our model should yield to 

a mathematical analysis eventually to give us such clue. 



IV. SUMMARY AND PROSPECTS 

In this work, an extended standard map which has 
spatial and temporal asymmetry is formulated. We have 
outlined its phase space properties, and have studied the 
unidirectional diffusion, which is thought to be relevant 
for molecular and biological motors. There exist several 
studies of great mathematical rigor on the standard map 
[]l2| , [l3| . Our extended standard map may hopefully be- 
come a subject of similar research. One immediate exten- 
sion of the current model is to introduce the randomness. 
Namely, one may consider s not as a fixed quantity, but a 
fluctuating number for each n. This is a simple way to in- 
clude the thermal fluctuation into the model. Lastly, the 
study of quantum version of the current model may be of 
great interest. This is particularly so, because the quan- 
tum version of standard map [|l^, |l^ is again a standard 
tool in the study of quantum chaotic diffusion. 
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